Complete Characterization of Quantum-Optical Processes 
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The technologies of quantum information and quantum control are rapidly improving, but full exploitation 
of their capabilities requires complete characterization and assessment of processes that occur within quantum 
devices. We present a method for characterizing, with arbitrarily high accuracy, any quantum optical process. 
Our protocol recovers complete knowledge of the process by studying, via homodyne tomography, its effect on 
a set of coherent states, i.e. classical fields produced by common laser sources. We demonstrate the capability 
of our protocol by evaluating and experimentally verifying the effect of a test process on squeezed vacuum. 
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PACS numbers: 

Construction of a complex machine requires precise charac- 
terization of each component's properties. In electronics, this 
information is obtained from network analyzers, which mea- 
sure circuit response to simple oscillatory inputs, and reveal 
the device transfer function. Optical quantum technologies, 
which can be used to build quantum computers 110, precise 
metrological systems |2], and unconditionally secure commu- 
nication yD, have similar characterization requirements. In 
this context, we are interested in the process associated with 
a quantum circuit component, i.e. in being able to predict the 
transformation that an arbitrary quantum state will undergo 
when subjected to the action of the component. 

A quantum process £ can be represented by a completely 
positive, trace-preserving linear map (superoperator) on the 
linear space IL(H) of all density matrices over Hilbert space 
H. It can be expressed as a rank-4 tensor that relates the ma- 
trix elements of the output £[p) and input p states in some 
basis: 
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where summation is from 1 to dim H. 

Characterization of a process (known as quantum process 
tomography, or QPT) means finding all components of the su- 
peroperator tensor. It can be implemented by determining the 
output state for each of the (dimH) 2 elements of a spanning 
set of JL(H). Such a direct approach to QPT |3] was exper- 
imentally realized on one qubit teleportation [5], the Hamil- 
tonian evolution of vibrational states of atoms in an optical 
lattice 0], on a two-qubit controlled-NOT gate 10,1] and Bell- 
state filter OJ. As an alternative, ancilla-assisted QPT exploits 
an isomorphism between processes and states flQ|| and has 
been used to characterize a controlled-NOT gate fllll and a 
general single qubit gate (HI [HI] ; see fl4ll for a comparative 
review of ancilla-assisted QPT. 

Existing QPT suffers from serious shortcomings, includ- 
ing either the requirement of an unwieldy set of input states 
for direct QPT or a high-dimensional entangled input state for 
ancilla-assisted QPT; these shortcomings deleteriously affect 
scalability and restrict accessible systems to very low dimen- 
sion. In optics, QPT has been applied to processes on one 



and two dual-rail qubits, with post-selection based on photon 
coincidences projecting the input and output states onto these 
qubit subspaces. This approach cannot provide complete in- 
formation about a state or a process because optical losses, 
imperfect sources, detector dark counts, and other imperfec- 
tions lead to departure from the qubit subspaces. Post-selected 
tomography can only estimate the fraction of such phenomena 
by comparing the coincidence rate and the photon production 
rate H. 

We introduce a scheme that enables complete characteri- 
zation of a general quantum-optical process. We use optical 
homodyne tomography followed by maximum likelihood re- 
construction to obtain full information on the process across 
all photon number sectors and also the coherence between sec- 
tors. The state reconstruction algorithm provides an efficient 
method for compensating losses in homodyne detection lH3l . 
As inputs, we use only coherent states that are readily avail- 
able from a laser source, so our method can be easily scaled 
up. 

We experimentally test our approach by characterizing a 
quantum process that consists of a simultaneous absorption 
and phase shift. The reconstructed superoperator allows us to 
predict, with a fidelity of over 99%, the effect of the process 
on a squeezed vacuum. 

Our method is based on the fact that any density matrix 
can be represented as a sum of coherent states' density ma- 
trices fill [nil . Although such a representation (the Glauber- 
Sudarshan P function) may be highly singular, it can be ar- 
bitrarily closely approximated with a regular P function. By 
measuring the process output for many coherent states and ex- 
ploiting the linearity, we can predict the process output for any 
arbitrary state. 

The Glauber- Sudarshan decomposition of a quantum state 
p is given by 



p = 2 J P p (a)\a)(a\d 2 a 



(2) 
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and where P p (a) is the state's P function, a is the coherent 
state with mean position and momentum observables (x,p) = 
( V2Re a, a/2Ihi a). We use the convention [x,p] = i and in- 
tegration is performed over the entire phase space. Therefore, 
if we know the effect \a)(a\ i— > g(a) = £(\a) (a\) of the pro- 
cess on all coherent states, we can predict its effect upon state 
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D p (a)g(a) d a. 



(3) 



An obstacle to direct application of this approach is posed 
by singular behavior of the function P p (a). Indeed, the P 
function exists only as a generalized function, more singular 
than the Dirac delta function, when the corresponding quan- 
tum state has nonclassical features iflill . 

This can be overcome by applying a theorem proven by 
Klauder II19I1 : for any bounded operator p there exists an op- 
erator pl with continuous and rapidly decreasing P function 
arbitrarily close to p in the trace-class norm. The Klauder 
approximation is obtained as follows: we assume that the 
Wigner function of p belongs to the Schwartz class <S 2 , i.e. is 
infinitely smooth and rapidly decreasing (which is the case for 
all physically plausible density matrices). The Fourier trans- 
form of the operator's Glauber- Sudarshan function P p (a) can 
be expressed as II 1 811 



Pp{kxi kp) — Wp(k X i kp) exp 




(4) 



where W p (k x , k p ) is the Fourier transform of the operator's 
Wigner function. The function defined by Eq. [4] always ex- 
ists and is infinitely smooth (albeit not necessarily square in- 
tegrate). We multiply P p (k XJ k p ) by a regularizing function 

G L (k x ,k ) = e -lf( k x - L )+f(- k x - L )+f( k P- L )+f(- k P- L )) 

(5) 

with f(y) = y 4 exp(-l/y 2 ) for y > 0, f(y) = for y < 0. 
This regularizing function is equal to 1 in a square domain 
of side 2L and rapidly drops to zero outside. The prod- 
uct PL, P (k x , k p ) = P p (k x , k p )G(k x , k p ) now belongs to the 
Schwartz class. Applying the inverse Fourier transform, we 
obtain the new Glauber- Sudarshan function PL yP (a), which 
defines the Klauder approximation fa. By choosing L suffi- 
ciently high [20], the operator fa can be made to approximate 
p arbitrarily well (Fig. 4A). 

As an example, we applied the Klauder approximation 
to squeezed vacuum, a nonclassical state characterized by a 
highly singular P function whose Fourier transform grows ex- 
ponentially with increasing k x and/or k p . We tested our pro- 
tocol with a state that has a noise reduction in the squeezed 
quadrature of —1.58 dB and excess noise in the orthogonal 
quadrature of 2.91 dB. The function P(k x , k p ) was calculated 
from the state's density matrix according to Eq. [4] and subse- 
quently regularized as described above using L = 5.2. Fig. 
1A shows PL(k X jkp) calculated from our experimental data 
and Fig. IB displays its inverse Fourier transform Pl((x). In 
Figs. 1C and ID we compare the Wigner functions of the orig- 
inal state and the one obtained from the regularized P func- 
tion. The two states exhibit a quantum fidelity of more than 
0.9999. 

Although the above method permits finding the process out- 
put for an arbitrary input state, it requires one first to deter- 
mine the input state's P function. This step can be avoided 
by calculating the process superoperator in the Fock basis, so 
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FIG. 1: Regularized Glauber-Sudarshan decomposition of the 
squeezed state. (A) Absolute value of the regularized Fourier trans- 
form of the squeezed vacuum P function. (B) Approximated P func- 
tion calculated from the inverse Fourier transform of PL,p(k Xl k p ). 
(C) and (D) Wigner representations of, respectively, the measured 
and approximated squeezed vacuum states. 



the output can be found from the input density matrix accord- 
ing to Eq. [B To this end, we express the Glauber-Sudarshan 
function as 



P p (a) 
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where P nrn (a) is the P function of the operator \n)(m\. 
We now replace these functions by their regularized versions 
PL,nm(ot) and rewrite Eq. [3]as 
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(a)g(a)d a, 



from which we determine the process superoperator as 
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£/T = 2 / P L ,nm(a)g lk (a)d z a. 
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Prior to applying the latter result to experiments, a num- 
ber of practical issues have to be addressed. First, parameter 
L must be chosen to ensure proper approximation of input 
states. The second issue is that in a realistic experiment, the 
measurement can be done only for coherent states whose am- 
plitude does not exceed a certain maximum a max . Finally, the 
experiment can only be performed with a finite, discrete set of 
coherent states. Density matrix elements Qik(&) for an arbi- 
trary a, required for calculating the superoperator, must then 
be obtained by polynomial interpolation. These matters are 
discussed in 12011 . 

A simplification arises for phase- symmetric processes, in 
which there is no phase coherence between the "processing 
unit" and input states. In this case, if two inputs p and p\ are 
different by an optical phase shift U(ip), the states S(p) and 
£(pi) will differ by the same phase shift: 



(9) 
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Then, if we know the effect of the process on a coherent state 
\a), we also know what happens to lae*^), so it is enough 
to perform measurements on input coherent states with real, 
positive amplitudes. For the process superoperator in the Fock 
basis, the phase symmetry implies that £™ n vanishes unless 
k — I = m — n. 

The process studied in our experiment was electro-optical 
amplitude and phase modulation of the optical field. The 
process was implemented using an electro optical modulator 
(EOM) followed by a polarizer. The field experienced mini- 
mal distortion when a bias voltage V\ = 100 V was applied to 
the EOM. Switching the voltage to V2 = 50 V produced bire- 
fringence, and thus losses at the polarizer, along with a phase 
shift. 

A continuous-wave Ti: Sapphire laser at 795 nm was the co- 
herent state source used for the device characterization. We 
reconstructed the input and output states at 1 1 different input 
amplitude levels between ai = and an = 10.9. In order 
to keep track of the relative phase shift, the EOM voltage was 
switched between V\ and V2 every 100 /is (Fig. 2A, top) while 
the phase of the local oscillator was linearly scanned by a 
piezoelectric transducer at 100 Hz. The homodyne photocur- 
rent was recorded with an oscilloscope. To obtain quadrature 
measurements, the photocurrent was integrated over time in- 
tervals of 20 ns. The bottom plot in Fig. 2 A shows the re- 
covered quadrature values after normalization to the vacuum 
noise. The time dependence of the local oscillator phase was 
recovered from the slow, sinusoidal variation of the average 
homodyne photocurrent as a function of time. 

In this manner, for each input amplitude, we sampled 
50,000 phase and quadrature values for both the input and 
output states and used them to calculate density matrices by 
likelihood maximization tl5l |2lj] (Fig. 2B). The output state 
reconstruction showed a phase shift of 36° and a loss of 34 % 
with respect to the input state. 

The interpolated experimental data have been used to de- 
termine the process superoperator tensor. We used the phase 
symmetry assumption in Eq. which is justified by the fact 
that the EOM driver is independent from the master laser. The 
elements £ ^ m of the tensor in the photon number basis are 
plotted in Fig. 3 A. This plot should be interpreted as follows: 
for a given input Fock state \m), the values of give the 
diagonal elements of the output density matrix. For example, 
the single-photon state 1 1) after passing through the EOM will 
be transformed into a statistical mixture of the single-photon 
and vacuum states. A theoretical prediction for the process 
tensor has been calculated using the Bernoulli transformation 
to account for a lossy channel and a phase shift superopera- 
tor; the superoperator diagonal elements in the Fock basis are 
displayed in Fig. 3B, and these diagonal elements bear close 
resemblance to the experimental result. A similar agreement 
was also obtained for non-diagonal terms of the superoperator, 
but it is not shown here. 

For additional verification, we applied this result to pre- 
dict the effect of the device on the squeezed vacuum state 
described in the previous section. This state was produced 
by pumping an optical parametric amplifier (OPA) in bow 
tie configuration with the second harmonic of the Ti: Sapphire 
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FIG. 2: (A) Time-dependent quadrature values acquired from homo- 
dyne detection of a coherent state with input a = 8.3. Black dots 
correspond to the state before the process; red dots, after the process. 
The top curve shows the EOM driving voltage. (B) The Wigner func- 
tion of the coherent state before and after the process. (C) and (D) 
Wigner representations of the measured output squeezed state com- 
pared to the one obtained from process tomography. 
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FIG. 3: The "diagonal" values of the superoperator £™ k m . (A): as 
obtained in the experiment. (B): theoretical model. 



laser and using a periodically poled KT1OPO4 crystal as non- 
linear mediurn|l22l [2311 . 

The state before (Fig. 1C) and after (Fig. 2C) the process 
was reconstructed using homodyne tomography as described 
above. By applying the process superoperator to the input 
squeezed state, we predict the process output (Fig. 2D). The 
maximum (minimum) quadrature noise variance amounted 
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to 2.19(-1.07) dB for the measured state, and 2.15(-0.95) 
dB for the prediction, corresponding to a quantum fidelity of 
0.9935±0.0002. 

Whereas here we demonstrate our tomographic method for 
single-mode inputs, multimode or multichannel processes can 
be characterized using multimode P representation, multiple 
homodyne detectors and feeding product coherent states as 
inputs. Theoretical and experimental analysis of the multi- 
mode setting will be discussed elsewhere. 

Our method overcomes significant limitations of previ- 
ous optical QPT schemes. Process characterization is not 
restricted to a Hilbert space associated with a specific qubit, 
and thus reveals the imperfections of a quantum information 
processing unit. Additionally, it uses only coherent states as 
inputs, which are readily available from the laser and whose 
intensities and phases are easily manipulated. This permits 
characterization of complex processes employed in quantum 
information processing and communication. 

Supporting Material. In order to compute the P function 
Pl, p {ol) by the Klauder approximation II 1 911 and the superop- 
erator terms from 
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(a)gik(a) d 2 a, 



(10) 



where Pi, jnm (a) is the regularized P function of the operator 
\n)(m\, the values of two parameters must be set: the size L 
of the Fourier domain, and the maximum coherent state am- 
plitude Qfmax- As evidenced by Fig. 4 A, no value of L can 
ensure universally high fidelity for all optical states. How- 
ever, for many practical applications it is reasonable to restrict 
the input states to subspace = span(|0) , . . . , \N)) with 
a limited number of photons. Under this restriction, we can 
seek the lowest L that, in the worst case over all density ma- 
trices in Hat, yields an approximation to the input state with 
required fidelity F. 

We employed the genetic optimization algorithm ll24f1 to 
calculate the required value of L for several low values of N 
and ascertained that the worst case scenario corresponds to the 
Fock state \N). This is not surprising, because the P function 
of this state contains the highest order derivative of the Dirac 
delta in Hat. We assume this rule to hold for arbitrary N. 

Restricting the amplitude of coherent states for which the 
measurements are performed effectively introduces finite in- 
tegration limits in Eq. 1, entailing additional fidelity loss. 
Fig. 4B shows the lowest a max required for approximating the 
TV-photon Fock state (which, again, appears to be the worst 
case in Hat) with a 99 per cent fidelity. Notably, this quantity 
is much larger than that expected from the behavior of the cor- 
responding Wigner function. This is because the regularized 
P function exhibits strong oscillations in a much wider phase 
space region than does the Wigner function. 

Even with a finite a max , reconstructing the superoperator 
requires knowing the process output for a continuum of coher- 
ent states. In an experiment, we choose a finite set of inputs 



\ot.i) from a laser source with varying, but known, amplitudes 
and phases, and perform homodyne tomography to determine 
the corresponding output density matrices gik (on) in the Fock 
basis. For an arbitrary a, density matrix elements Qik(ot) are 
then obtained by polynomial interpolation. 

The degree of the polynomial over a can be reduced for 
Gaussian and near-Gaussian maps, which comprise all linear- 
optical processes and nonlinearities. For each measured 
output Q(a.i) we infer its mean complex amplitude (a) = 
Tr (dg(ai)). Then we apply the displacement operator 



(a) 



D(-(a))g(a i )D^(-(d)). 



(11) 



to center the Wigner function at the origin of the phase 
space. The matrix elements g$ , as well as the mean am- 
plitude (a), can now be fitted with a low-degree polynomial 
over a, allowing for a highly efficient reconstruction process. 
Good performance of this algorithm for Gaussian preserving 
maps is explained by the fact that Gaussian states are en- 
tirely characterized by the first and second centered moments 

(x - y/2 Re(a)) m (p - V21m(a)) n of their Wigner distribu- 
tion. 
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FIG. 4: Requirements for faithful approximation of photon number 
states with a regularized, restricted Glauber-Sudarshan decomposi- 
tion. (A): the low pass filtering parameter L needed to achieve a 
99.99-% fidelity (•) and 99-% fidelity (o). (B): the lowest a m ax re- 
quired for a 99-% fidelity. 
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